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Abstract 



Using periodic orbit theory, we evaluate the form factor of a quantum graph to which a 
very weak magnetic field is applied. The first correction to the diagonal approximation 
describing the transition between the universality classes is shown to be in agreement 
with Pandey and Mehta's formula of parametric random matrix theory. 
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1 Introduction 



Since more than a decade, it has been known that the energy level statistics of classi- 
cally chaotic systems universally follows the prediction of random matrix theory. As the 
semiclassical theory of chaotic systems is described in terms of classical periodic orbits, 
it is preferable if we can also understand the universal level statistics in terms of them. 
However, until recently, such an understanding had been limited within the framework of 
Berry's diagonal approximation for double sums over periodic orbits]!]]. 

A recent breakthrough was brought about by Sieber and Richterp], §]. They identified 
the pairs of periodic orbits giving the first off-diagonal correction and showed that such 
pairs indeed exist in chaotic systems. The first off-diagonal correction was then found to 
be in agreement with random matrix theory. Berkolaiko, Schanz and Whitney rederived 
the same correction term for quantum graphs [|]] and further succeeded in evaluating the 
second correction ||. 

In this paper, we investigate a quantum graph and extend the calculation of the first 
off-diagonal correction to the case intermediate between the absence and presence of a 
magnetic field. In random matrix theory, this intermediate case corresponds to Pandey 
and Mehta's two matrix model|| |7|]. Chaotic systems with time reversal symmetry (with- 
out magnetic field) are described by GOE(Gaussian Orthogonal Ensemble) of random 
matrices. On the other hand, when time reversal symmetry is broken by the application 
of a magnetic field, GUE(Gaussian Unitary Ensemble) becomes suitable. Pandey and 
Mehta formulated an intermediate random matrix ensemble (two matrix model) between 
GOE and GUE and were able to derive the correlation functions among the eigenvalues. 
Within the diagonal approximation, Bohigas et al.|| already showed that the periodic 
orbit theory was in accordance with the intermediate random matrix ensemble. We will 
show that, by applying Berkolaiko, Schanz and Whitney's method, the first off-diagonal 
correction is also in agreement with Pandey and Mehta's formula. 

Let us explain the model. We consider quantum mechanics on a graph |9|, [10], [IT], |12|j . 
A graph consists of vertices connected by bonds. For example, in a globally coupled 
graph, every vertex is connected to all vertices, while in a star graph, one particular 
central vertex is connected to all of the others. 

On the bonds a particle behaves like a free particle and on the vertices it is scattered 
according to a given scattering matrix. On the bond (j, I) connecting the j-th vertex and 
Z-th vertex, the Schroodinger equation 



holds, where Aji is a magnetic vector potential satisfying Aji = —Ay. Therefore, on the 
bond (j, I), the wave function ^(xji) is proportional to 



The relative amplitude of this wave function is decided by the scattering matrices on the 
vertices. Let us suppose that, if a wave function has the amplitude 1 on the bond (j, I), 




(1.1) 



exp (ikxji + iAjiXji) . 



(1.2) 



then the amplitude of the wave function on the bond (/, m) is cr^ m . Then the scattering 
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matrix connecting the bond (J, I') and (l,Tn) is 



Sji',im = af^exp^kLim + iAimLir^di^, (1.3) 

where Lji is the length of the bond (j, I) satisfying Lji = L\j. The total number of the 
vertices and directed bonds are denoted by iV and B, respectively. In this paper, we 
consider only a globally coupled graph (every vertex is connected to all vertices), so that 
B is equal to N 2 . We further assume that the scattering matrix is given by the discrete 
Fourier transform 

*S, = -^e 2 ™^. (1.4) 



Chaotic dynamics in the classical limit is characterized by an approach to an equidistri- 
bution over all bonds. Noting that an analogue of the classical Frobenius- Perron operator 
is given by M m /^ m = \S m 'i,i m \ 2 , we can write it as 

lim M t m , l lm = l/B. (1.5) 

t — >oo 1 



For the scattering matrix ( |1.4j) , ( |1.5| ) holds even for finite t. Berkolaiko, Schanz and 



Whitney]^] argued that, if the convergence of ( |1.5| ) in the limit t — > oo is sufficiently fast, 
the form factor is in agreement with random matrix theory. We expect that a similar 
condition holds in the intermediate case. However it is not discussed here and left for 
future works. 



2 Periodic Orbit Theory 

Let us define P and Q as sequences of vertices [pi,P2, ■ ■ • ,Pt] and [q%, q 2 , ■ • • , q t ], respec- 
tively. Each of pj and qj takes the values 1, 2, ■ • • , N. According to Berkolaiko, Schanz 
and Whitney H, the form factor K(t) can be written in terms of the scattering matrices 



as 



k(t) = 4(l Tr5 T) 



where 



and 



B 
1 

B 



- lim k^Y /^fca^ ) e iL w^(fc+^iP2) a (P2) e ^2P3( fc +^ 
B K - > °° p Q Jo PtT2 P1P3 



i) 



i<T (Pt) JLp tP1 (k+Ap tP1 ) ( qi )* -iL nq2 (k+A qiq , 2 ) _ _ _ (q t ) * „-iL qtqi {k+A qtn ) 
u pt-ipi c u qtq2 c U Qt-iqi c 



_L \ A „ A * p'^(^ j piP2 j ^piP2^ ^Lp tP1 A PtP1 L qiq2 A qiq2 ■■■ L qtqi A qtqi \ r 



P,Q 



(2.1) 



= ^•••ofcU L P = L piP2 + ... + L PtP1 , (2.2) 
= ^ 1 ql--^i qv L Q = L qiq2 + ... + L qtqi (2.3) 

r = t/B. (2.4) 
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We are interested in the scaling limit t — > oo, B — > oo with r fixed. Then we expect that 
the diagonal term and its corrections give the expansion around r = 0. 

The diagonal term comes from the cases in which circular permutations of P and Q 
(or the reversal of Q) coincide. It is explicitly written as 



K diag (r) = -Lj2\^p\ 2 [e 2i ^ Lpm+lApjPj+1 +^\, (2-5) 
p 

where p t +\ = Pi- Ori the other hand, the first off-diagonal correction K (r) comes from 
the pairs of self-intersecting orbits differing in the orientation of a single loop. We suppose 
that the orbits are self- intersecting at a vertex a = p\ — p t >. Then, in the first off-diagonal 
correction, the orbit sum is taken over the pair P = [a, l±, a, l 2 ] and Q = [a, li, a, l 2 ]. Here 
h = [P2,P3, • • • ,Pt'-i], h = \p t '+i,Pt'+2, ■ ■ -,Pt], respectively, and l 2 = \p t ,Pt-i, • • • ,Pt'+i] is 
the reversal of the sequence l 2 . That is, K oS {r) can be written as 

.2 t-2 



K otl (r) 



-° t'=4 P 



(a) (c*Wa)* (a)« 





2 . 


..UW pUW 12.. 


■W (d) 1 






1 p t ,_ 2 a\ \ w ap t , +2 \ 




exp 


2i 


^Lda-^da ^Pt-iPt^Pt- 


-lpt + ■ ' ■ 



(2.6) 



Here a, b, c, d are identified with p2,pt>-i,Pt'+i,Pt, respectively. As indicated by the prime, 
the sum over P obeys the restriction that l\ is not identical to the reversal of itself (h ^h). 

We are now in a position to calculate the diagonal contribution and the first off- 
diagonal correction. For simplicity, we set Lji = 1/2 from now on. 

(1) the diagonal contribution 



— \ A p\ 2 [f: 2 '^ 7 '"'"'- 1 " 1 "'"'-' + 1 
B 



p 

L v — 

BpN 1 



o 2l Hj L PjPj + l A Pj P] + l _|_ I 



= LJ_\ N t + Tr t 

B N l L y 

Here the elements gji of an N x N transfer matrix g is defined as 



9ji 



1, 3 = I, 



(2.7) 



(2.8) 



Since we are interested in the long time limit t — > oo, we only need to know the behavior 
of the largest eigenvalue A of g. For small Aji, the largest eigenvalue A can be evaluated 
by the perturbation method. Let us decompose g as 



9 = 9o + 9 , 



(2.9) 



where go is an A" x N matrix with all the elements equal to 1. Then the largest eigenvalue 
of g is A" and the corresponding normalized eigenvector x x is (l/\/N, 1/y/N, ■ ■ ■ , l/\/N) T 
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(Here a superscript T means a transpose). All the other eigenvalues of go are and let us 
denote the corresponding orthonormalized eigenvectors by x 2 , x 3 , • • • , x^v- Then we regard 
g' as a small perturbation and apply the standard technique of perturbation theory to 
evaluate the largest eigenvalue. The second order perturbation gives 



1 N 



n=2 

which yields 



A ~ N + xf </xi + ^ E l x n Ail 2 , (2.10) 



A-iv-^E4 + 4El x n« x il 2 ' (2- 11 ) 

iV j<l iV n=2 



where a is an N x N matrix with the elements 



a,-,= { °' J J? (2.12) 



In the limit N — > oo, the sum C of the contributions from the second and third terms 
is normally O(N) (for fixed Ajj). We then generally write C ~ — with b proportional 
to the magnetic field applied to the graph. For example, if J2i^j Aji = for all j and 
\Aji\ = A for all j, I (this is possible when N is odd), the third term vanishes and the 

second term gives C = —{N - I) A 2 / 2 NA 2 /2, so that b = A 2 / 2. The trace of the 

multiples of the matrix g can be estimated as 

Trg* ~ A* ~ (iV — Nbf ~ iV'e -6 * (2.13) 

with t — » oo, 6^0 and 6£ fixed. Putting this, we arrive at 

^ diag (r) ~ 4( 1 + e ~ 6 *)- (2-14) 
B 

This is the diagonal contribution to the form factor in a weak magnetic field, as we 
shall see below, in agreement with random matrix theory. We now proceed to the first 
off-diagonal correction. 

(2) the first off-diagonal correction 

As noted before, the contribution K°^(t) from self-retracing orbits with l\ = l x is removed 
from the periodic orbit sum (|2.6|) for K oS (r), since it is already included in the diagonal 
contribution. In spite of that, it is convenient to first consider the sum with no such 
restriction. The periodic sum including the self-retracing orbits can be readily evaluated 

as 

K oS (r)+K^(r) = t g AT*'-*- £ (1 _ 5cd ) g ac 

t'=4 aabcd 

= | eV-- 2 EMi -Us, 

t'=4 acd 

Here we utilized the unitarity of the matrices <j( a \ As to the contribution from the self- 
retracing orbits, if we neglect completely self-retracing paths which are exponentially few, 



/ \ cd gda<J da °bc a ca °bd 



J-t'-l 



c,d 



g da = 0. (2.15) 
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a diagrammatic cancellation leaves only the terms with t' = 4 and t' = 5, in which the 
factors 1 — 5 c d are removed 01. Therefore we find 



K°\r) = -K°°(t) 



-off/ 



B 



aacd 



J-5 



c,d 



9do 



aacd 



c,d 



9do 



B 



iV 1_ *Tr ( g l ~ 6 ) + iV^Tr ( g 



t-3 



2-tm_ t-A 



(2.16) 



Substitution of ( |2.13|) into ( |2 . 1 6| ) leads to 



K oS (t) ~ 



-bt 



[2.17) 



As a result, the diagonal and the first off-diagonal terms are summed up to yield 



K dia Hr) + K oS (t) ~ r + re" 6 ' - 2rV w 



(2.18) 



3 Random Matrix Result: Pandey-Mehta Formula 

In this section we calculate the Fourier transform of random matrix result (Pandey-Mehta 
formula) to derive a prediction of the form factor K{r). Let us define the Fourier transform 
of the two energy level correlation function Y(r; p) as Y(k; p). The form factor K(k) can 
be then written as 

K(k) = l-Y(k;p). (3.1) 

Pandey and Mehta's formula || is 

Y{r-p)= sinTrrX 2 _ 1 r*^ r°° ^ / M ^ sin ( fc2r ) e 2p 2 (^ 2 )(^ 2 ). (3 . 2 ) 

Here r is the distance between the two energy levels and p is a parameter corresponding 
to a weak magnetic field. In the limit p — > 0, V(r; p) becomes the two level correlation 
function of GOE, while, in the limit p — > oo, y(r; p) approaches that of GUE. 

The Fourier transform Y(k; p) can be evaluated as 

Y(k; p) = 1 - k — r dkt-^e- 2 " 2 ^^, 0<k<n, (3.3) 

2n Jn-k ki + k 

where k = 2nk. Let us define 

F(k) = J- f djfei-A^e-^fci+K) (3 .4) 

27T Jvr-fc fcl + k 

with c = p 2 /c fixed. Then, for small k, F(k) can be expanded as 

F(k) = F(0) + F'(0)k + V'(0)P + • • ■ , (3.5) 
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where F'(k) and F"(k) are the first and second derivatives of F(k), respectively. The 
derivatives are readily evaluated as 

no) = ^e" 4 - 

F"(0) = -^e" 4 ™, (3.6) 

which lead to (for small k) 

Y(k-p) ~ l — k — F'(Q)k — ]-F"(0)k 2 

= l — k — ke~ ilTC + 2/c 2 e" 4?rc . (3.7) 

Therefore we arrive at 

tf(Jfe) ~ k + ke~ 4nc - 2k 2 e~ 4nc , (3.8) 
which is in agreement with periodic orbit theory with an identification 

t = k, bt = 4vrc = 8n 2 p 2 k. (3.9) 

In summary, using the periodic orbit theory we evaluated the first off-diagonal correc- 
tion to the form factor of a quantum graph in a very weak magnetic field. It was found 
that the result was consistently in agreement with Pandey-Mehta formula of random ma- 
trix theory. The extension to the second off-diagonal correction seems promising, since it 
was already worked out in the absence of a magnetic field || . 
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